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We show that the novel features observed in Raman 
experiments on optimally doped and underdoped Bi-2212 
compounds in B\ g geometry can be explained by a strong 
fermionic self-energy due to the interaction with spin fluctua- 
tions. We compute the Raman intensity R(u) both above and 
below T c , and show that in both cases R(uf) progressively de- 
viates, with decreasing doping, from that in a Fermi-gas due 
to increasing contribution from the fermionic self-energy. We 
also show that the final state interaction increases with de- 
creasing doping and gradually transforms the 2A peak in the 
superconducting state into a pseudo resonance mode below 
2A. We argue that these results agree well with the experi- 
mental data for Bi-2212. 

keywords: A Superconductors, E Inelastic Light Scattering, 
D Electron Electron Interactions. 

The form of the fermionic spectral function in opti- 
mally doped and underdoped cuprate superconductors 
has been the subject of intensive experimental and theo- 
retical studies over the last few years. ARPES and neu- 
tron experiments demonstrated that the fermionic spec- 
tral function undergoes a substantial evolution with de- 
creasing doping, and for underdoped cuprates is very dif- 
ferent from the one in a Fermi-gas (FG) Q. Raman scat- 
tering in Big geometry is another spectroscopy to study 
the electronic properties in this momentum region. For 
electronic Raman scattering, the intensity R(u>) is in gen- 
eral given by the imaginary part of the fermionic particle- 
hole bubble at small external momentum q and finite 
frequency, weighted with the Raman vertices ^-^]. In 
this paper we focus on the Bi g Raman scattering. In 
Big geometry, the Raman vertex Vb 19 oc cosfc x — cosfc y , 
and the scattering thus mostly probes the vicinity of 
(0, 7r) ||. Recent Raman experiments on overdoped, op- 
timally doped and underdoped Bi— 2212 demonstrated 
that the Bi g Raman intensity R(uj) undergoes significant 
changes with decreasing doping and progressively 

deviates from predictions of a FG theory. Indeed, accord- 
ing to the FG theory, R(u>) in the normal state is finite 
only for nonzero external momentum q, and vanishes at 
to > vpq, where vp is the Fermi velocity flo|| . In the 
superconducting state, R(u>) is finite due to a particle- 
hole mixing and possesses a peak at oj — 2A where A 
is the maximum of the superconducting gap A^ ||. Be- 
low 2 A, the intensity scales as cu 3 at small frequencies 
due to the presence of nodes in the d— wave gap ||. In 
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FIG. 1. The Bi g electronic Raman scattering spectra at 
T = 4 K for overdoped (OD), optimally doped (OP) and un- 
derdoped (UD) Bi 2 Sr 2 CaCu208± 4 (Ref. [§). The inset shows 
the high energy continuum in the normal state for the under- 
doped T c = 60 K sample and its fit to the a; 1 / 2 behavior. 



Fig. [j] we present the experimental Raman intensity for 
several Bi— 2212 compounds both below and above T c 
||. In the overdoped, T c = 82K material, the behav- 
ior of R(oj) is qualitatively consistent with the FG the- 
ory, i.e., it is featureless above T c (not shown), while at 
T < T c , R(uj) tx lo 3 at small frequencies JOJ], and exhibits 
a sharp peak at 2A « 50 meV. At and below optimal dop- 
ing, however, the form of R{lo) is inconsistent with the 
FG theory. In the normal state, the intensity increases 
with frequency as R(lu) ~ -y/uJ, and saturates at a few 
hundred meV [f^Ji"l||. (see the inset in Fig. |l|). In the 
superconducting state, the key experimental observation 
is that with underdoping, the peak in R(u>) occurs at 
progressively smaller frequency than the " 2 A" extracted 
from tunneling experiments |]l2|| , and almost saturates at 
about 75 meV Q . Simultaneously, the low frequency be- 
havior becomes predominantly linear in to with a kink 
around 40 — 50 meV, while above the peak the inten- 
sity develops a dip at about 90 meV and at even larger 
frequencies recovers the normal state yfuj form. 

In this Letter, we make two points. First, we argue 
that the inconsistency of the experimental results for 
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FIG. 2. The full B\ g Raman intensity in a superconductor 
(a), and the full B\ g Raman vertex (b). The wavy line and 
dashed rectangle are the fully renormalized effective interac- 
tion in d— wave and s— wave channels, respectively. 



Bi-2212 at and below optimal doping with the FG sce- 
nario indicates that the fcrmionic self-energy is large and 
substantially modifies the form of the fermionic Green's 
function both in the normal and in the superconduct- 
ing state. Phenomenologically, this effect was considered 
in fl3[| . Here we perform a microscopic analysis within 
a spin-fcrmion model in which the fermionic self-energy 
emerges due to an interaction with overdamped spin fluc- 
tuations and is the largest for fermions in the vicinity 
of (0, 7r). Since the Raman intensity in B\ g geometry 
predominantly probes the region near (0,7r), it directly 
reflects the changes in the fermionic spectrum imposed 
by the large self-energy. Second, we show that a mag- 
netically induced final state interaction replaces the 2A 
peak in R{u>) by a pseudo resonance peak at a smaller 
frequency tv res . Near optimal doping, this effect is al- 
most unobservable, but it becomes visible in underdoped 
cuprates and explains the discrepancy between the Ra- 
man peak frequency and the one in the density of states. 

We now turn to the calculations. The Raman intensity 
in a superconductor is given by a set of fermionic bubbles 
made of normal and anomalous Green's functions [ fbtUlB] 
(Fig.||). We first compute R(w) without final state inter- 
action. In this approximation, only the first two diagrams 
in Fig. ||a contribute (with bare vertices), and we have 

R{u) txlm J dk dfl Vf lg (k) 

x(G sc (k,n+) G sc (k,n-) + F(k,n + ) F(k,Q-)), (i) 

where fl± = £1 ± cj/2 and 

G sc (k, oj) = G-\-k, ~ W )/(G- 1 (fc, lu) G~ l {-k, -u) + A 2 k ) 
F(k,u) = iA k /(G- 1 (k,u) G-^-^-^+Al) . (2) 

Here A& is the d-wave superconducting gap, and 
G~ 1 {k,tu) = lu — efe+5 2 S(k, lu), where g is a dimensionless 
spin-fermion coupling, and S(k,w) is the A-dependent 
fcrmionic self-energy. Theoretically, g 2 oc where 
£ is the magnetic correlation length p6| . In strongly 
overdoped cuprates, g < 1, and the system behavior 



resembles that in a FG. However, at and below opti- 
mal doping, g > 1 in which case the self-energy over- 
shadows the bare ui term in the Green's function, i.e., 
^fcw) n g*X(k,w) -e k . 

The form of the fcrmionic self-energy is an input for 
our Raman calculations. Obviously, the self-energy is the 
largest near (0, tt) and symmetry related points where 
the scattering by nearly antiferromagnetic spin fluctua- 
tions is the strongest. The self-energy in this k— range 
also strongly depends on doping as evidenced by ARPES 
data. Since our goal is to relate the doping dependent 
changes in R(lu) with those in E(fc,cj), we restrict our 
consideration to the vicinity of (0, tt) where A k is close 
to its maximum value A. We will argue below that for 
strong coupling, the B lg Raman intensity is dominated 
by k near (0, tt) down to lu <C A, and crosses over to 
lu 3 behavior due to the nodes of the d-wave gap only at 
vanishingly small frequencies. 

It has been argued ||l6| , |l7f that for 9 > 1, Y,{k,Lu) 
near (0, tt) is independent of the quasiparticle energy 
up to corrections 0((\ogg)/g 2 ), and in a superconduc- 
tor behaves as T,(lu) oc u> at lu <C 2A, and as oc 
e l7T ^ 4 ^/\uj\sgnuj at w > max(2A, lu s j), where A = A/g 
is a measured gap, and Lu s f oc £~ 2 is a typical relaxation 
frequency of spin fluctuations. The physical reasoning 
here is twofold. First, in the normal state, the scat- 
tering by nearly-critical overdamped spin fluctuations 
yields S(w) = 2u/(l + ^/l — i\uj\/uj s f) which displays 
a crossover from a Fermi-liquid behavior at to < tu s f to 
a quantum-critical behavior for to > uj s f where S(cj) oc 
\fui Second, in a superconductor, the strength of the 
scattering by spin fluctuations is reduced below 2A due 
to a feedback effect on the spin damping jlTj . The calcu- 
lation of the full is rather involved and requires one 
to solve a set of two coupled complex integral equations 
for the fermionic self-energy and the spin polarization 
operator. Below we use the approximate self-consistent 
solution of this set of equations |17]] which correctly repro- 
duces the behavior of E(cj) at large and small frequencies. 
In the latter case it yields ~ ui, i.e., the peak in the 
spectral function occurs right at lu — A. 

We now proceed with R(ui). We assume that the 
density of states near (0, tt) depends only weakly on e k 
and replace the k— integration in (Q) by the integration 
over ek- We then obtain (neglecting overall prefactor) 
R(lu) = Imx{oj) where 
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(3) 



and D{n±) = y/X 2 {n ± ) - A 2 . 

The constant C > is a real number. Its presence in 
(^) is related to the fact that Eq. (|l|) with fc— integration 
substituted by the integration over lacks convergence. 
In this case, one cannot simply interchange frequency and 
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energy integrations and has to include a regularization 
procedure 18 1. The value of C is irrelevant for the calcu- 
lations of the Raman intensity without vertex corrections 
as it does not contribute to Jmx(a>). It however becomes 
relevant when one includes the effects of the final state 
interaction which accounts for the renormalization of the 
Raman vertex (see below). These vertex corrections in- 
volve the d— wave component of the effective interaction, 

which decreases at large frequencies and therefore pro- 
vides the physical regularization of the Raman bubble. 
In the spin-fluctuation approach, Yd starts decreasing at 
w ~ €p/g 0] ■ Another regularization is provided by the 
fact that the integral over is cut at ep . Obviously, the 
value of C depends on the ratio cf/S- For cf ^ g which 
is the case at weak coupling, the cutoff in a momentum 
space dominates. In this situation, an adequate way to 
evaluate the Raman bubble is to integrate first over fre- 
quency. This is how previous calculations of R(lu) have 
been performed ||[|. One then obtains (7=1. However, 
in the strong coupling limit ep "C <?, which is likely to be 
satisfied in underdoped cuprates, a more adequate way 
to evaluate the Raman bubble is to integrate first over 
as we did. In this case, (7 = 0. 

We now analyze Eq. (^) first in the normal state and 
then in the superconducting state. In the normal state, 
a substitution of the explicit form of into (j^) yields 
R(lu) = R(u„f) ^(uj/uj s }) where <fr(ai) » 1.07a; for x -C 1, 
and <&(x) « l.lZy/x for 1 < x < g A . We see that at 
lu > u,/, relevant to experiments, R(u) cx y/u). We fitted 
the normal state data by this form and found almost 
perfect agreement with the experiment (see the inset of 
Fig. 0). For even larger x > g 4 , the bare u> term in 
the quasiparticle Green's function begins to dominate, 
and the theoretical R(co) saturates, passes through a very 
broad maximum at x = 2.24g 4 and then slowly decays as 
1/y/x. This also agrees with the data. A similar behavior 
was found numerically in Ref. [|Ts|] . 

In the superconducting state, the form of R(w) clearly 
depends on the ratio b = A/ui s f. Experimentally, b < 1 
in strongly overdoped cuprates, and b 3> 1 in strongly 
underdoped cuprates |Q. For small b, Eq.(||) expect- 
edly reproduces the FG result for the Raman intensity: 
R(lo) = for oj < 2A and R(uj) cx (wVw 2 - 4A 2 )" 1 for 
u> > 2A. The inclusion of the k— dependence of the gap 
corrects this behavior at small frequencies where it yields 
a finite R{to) cx (ui/A) 3 M, and also very near 2 A where 
it changes the square root singularity to a logarithmic 
one. For large 6, however, the form of the Raman in- 
tensity substantially deviates from the FG result. We 
computed numerically the Raman intensity for different 
values of b, and present the results in Fig. || We em- 
phasize four key features in R(u>): (i) the 2 A peak is 
still present even for large b, (ii) with increasing 6, the 
peak frequency becomes progressively larger than 2A, 
(iii) there is a dip in R(u>) above the peak, (iv) at small 
frequencies, u> <§; 2A, R(u>) is predominantly linear in u. 
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FIG. 3. A calculated Raman intensity, Eqs. (^,^|) for dif- 
ferent values of b = A/uj s f. Larger b correspond to stronger 
self-energy corrections which destroy fermionic coherence. 
Observe the broadening of the peak with increasing b and 
the development of the dip at around 3 A. 



These features in the Raman intensity can all be under- 
stood analytically by analyzing Eq. (||) in various limits. 

The broadening of the peak, the dip and the linear 
behavior at low frequencies all agree with the features 
present in the data in Fig. [l]. However, one key dis- 
agreement with the data remains. Namely, we calculated 
the density of states N(oj) = j dkImG sc (k,uj) along the 
same lines as the Raman intensity and found that the 
peak in N(lu) is located at almost exactly a half of the 
peak frequency in R(u>), although both peak frequencies 
shift to larger values with increasing b. The data, we 
remind, show that with underdoping, the peak in R(u>) 
occurs at progressively lower frequencies than twice the 
peak frequency in N(u>). 

We now argue that the experimentally observed down- 
turn renormalization of the peak in R{lo) compared to 
twice that in N(ui) is due to a fact that the final state 
interaction between scattered quasiparticles gives rise to 
a bound state below 2A. 

The final state interaction in both s— wave and d-wave 
superconductors have been considered several times in 
the literature j3], |l^ , |l5| , pof . It gives rise to corrections to 
the particle-hole vertex (Fig. |2p), and also introduces an 
additional scattering process which mixes particle-hole 
and particle-particle channels (Fig. ||a). A sim- 

ple experimentation shows that the vertex corrections to 
the Big vertex involve the d— wave component of the ef- 
fective interaction in the zero-sound channel, while the 
"mixed" diagram (the last one in Fig. contains the 
fully renormalized s— wave vertex in the particle-particle 
channel. 

Several authors argued ]2(J that for a d— wave super- 
conductor with spin-independent interaction, the effec- 
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tive d— wave coupling in the zero-sound channel is re- 
pulsive, and the bound state does not appear unless the 
mixed diagram prevails over the conventional, RPA-type 
vertex renormalization. We, however, demonstrate below 
that for magnetically-mediated d— wave superconductiv- 
ity there is an additional sign change between vertices in 
the zero-sound and the Cooper channels, and this eventu- 
ally gives rise to attractive d— wave coupling in the zero- 
sound channel. 

Indeed, consider the effective interaction between 
fermions mediated by the exchange of spin fluctuations 
with momenta near Q. We have T = —g 2 x(q, fi) a a pa 1 f l 
where q = k — k' and Q = uj — uj' are transferred momen- 
tum and frequency, respectively. To simplify the discus- 
sion on the sign of the interaction, assume that T has a 
dominant d— wave partial amplitude Td, i.e., 



T(k - k',u - uj') w d k d k > T d (Japans 



(4) 



where dk are d— wave eigenfunctions. In general, Td 
depends on the transferred frequency. However, one 
can straightforwardly demonstrate that for a relaxational 
form of the spin susceptibility, the frequency dependence 
of Td becomes relevant at frequencies u> ~ ^> A [[f6| . 
The vertex renormalization on the other hand is mostly 
determined by much smaller uj ~ A. In this situation, Td 
can, to a good accuracy, be approximated by a constant. 

Let us first verify that there is a superconducting insta- 
bility in the spin-singlet particle-particle channel. Substi- 
tuting (||) into the equation for the full particle-particle 
vertex and making use of the identity 



0~a/30~j6 = T — 3S 



(5) 



where T,S — (SapS-^s ± S a sS 1 p)/2 are triplet and sin- 
glet spin configurations, respectively, we obtain (using 
for simplicity Fermi gas Green's functions) 



j^tot 

1 s 



-3r d /(i + 3r d l) 



(6) 



where L = log (uj max /T) > 0. We see that the behav- 
ior of the total particle-particle vertex in the spin singlet 
channel depends on the sign of Td- Several authors have 
demonstrated that near the antiferromagnetic instability, 
Td < [p2| . This obviously implies the superconducting 
instability. However, if the interaction were spin indepen- 
dent, then the d— wave instability would require a positive 
T d . 

Let us turn to the B\ g Raman intensity. The sim- 
plest way to check the sign of the vertex correction is to 
consider a Fermi gas in the normal state, and compute 
the density-density correlator at zero frequency and at a 
finite momentum. In the normal state, the mixed dia- 
gram does not contribute, and the ladder series of vertex 
correction diagrams can be easily summed up. We found 



i + T d Xd(q) 



where Vb 19 is the bare B\ g vertex, and 

Xd(q) = / dk (d k ) 



(8) 



is the uniform d— wave susceptibility (e± = ek±q/2, and 
Q(x) = 1 when x > and Q(x) = when x < 0). Ob- 
viously, Xd(q) > 0. We see that for negative Td, the 
Big vertex in the particle- hole channel is enhanced, i.e., 
the effective interaction in this channel is attractive. We 
show below that in the superconducting state this attrac- 
tion gives rise to a pseudo-resonance below 2A. Previous 
analytical studies [^0| obtained the same expression as in 
(||), but they considered spin-independent d— wave inter- 
action, and therefore set Td > 0. In this situation, vertex 
corrections reduce the B\ g vertex and do not give rise to 
a resonance behavior. 

The conclusion that vertex corrections reduce the B\ g 
vertex has recently been reached in numerical studies of 
the Hubbard model (|lj . The reasons for the discrepancy 
with our analysis are not clear to us because the effective 
interaction selected in J2l| is apparently also mediated 
by spin fluctuations. 

We now continue our analysis of the superconducting 
state. The vertex renormalization is given by a set of 
diagrams in Fig. |2|. For a spin-mediated interaction, the 
s— wave coupling is repulsive such that the last diagram 
in Fig. |^a does not contain a low-energy resonance mode 
and can be safely neglected. We are then left with the 
ladder series of vertex correction diagrams. For a fre- 
quency independent Td, the series of vertex corrections 
is then geometrical and yields 



Rfull ( w ) 



Imx(to) 



(9) 



(7) 



(l + T d Rex(uj)) 2 + (T d Imx(uj)) 2 

where, we remind, Td < 0, and x( w ) i s given by (|^). 
Now recall that in a superconductor Jmx(w) is small at 
lu < 2A. Evaluating Rex(ui) we found that it is posi- 
tive below u) < 2A. Thus for small frequencies, we found 
Re X (uj) = A(lu/A) 2 where in a FG, A = 1/3. Substi- 
tuting this result into (||), we find that Rf u ii possesses 
a resonance peak below 2A, at a frequency u) = ui res 
where \Td\ReR(uj res ) = 1. Actually, the solution for 
cores exists already at weak coupling simply because in 
a FG, Rex(u) diverges logarithmically as to approaches 
2A from below. As the coupling increases, the peak posi- 
tion progressively deviates downwards from 2 A. A simi- 
lar reasoning has been previously applied to explain the 
presence of the resonance peak in the neutron scattering 
data below T c @. 

Indeed, the bound state which we found is only a 
pseudo resonance because in a d— wave superconductor, 
Imx(u) is finite for all u> ^ because of the nodes of the 
gap. On the other hand, the very existence of the peak 
only requires a reduction of Imx(ui) below 2A which is 
a natural consequence of a reduction of the fermionic 
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FIG. 4. The fits of the theoretical Raman intensity with 
final state interaction, Rf u u(uj), Eq. [] to the experimental 
data for slightly overdoped, optimally doped and underdoped 
Bi2212 materials. We use the theoretical fact that in the 
absence of the final state interaction, the peak frequency in 
the Raman intensity, even at strong coupling, is almost ex- 
actly twice the peak frequency in the density of states, and 
extract an effective A from the photoemission and tunnel- 
ing data ||l|[l2|. The parameter Td measures the strength 
of the vertex corrections which, we argue, give rise to a 
pseudo-resonance in Rf u u(u). Observe that with underdop- 
ing, the peak in Rf u u(u}) progressively deviates down from 
2A. The presented fits are for b = 20. Fits using smaller b 
reproduce the shape of R(u>) but yield sharper peaks than in 
the data. 

spectral weight at low frequencies in a superconductor. 
Moreover, as the spectral weight reduction occurs already 
in the pseudogap regime |jj , our theory predicts that the 
Raman peak should survive above T c and disappear only 
at a temperature where the pseudogap behavior becomes 
invisible. 

In Fig. | we fitted the data from Fig. |l| to Eq. (|) 
using as an adjustable parameter. We found that in 
overdoped cuprates, the pseudo resonance frequency is 
almost indistinguishable from twice the peak frequency 
extracted from the tunneling data However, as the 
system moves towards lower doping, | | obviously in- 
creases, and the peak in Rf u u(uj) moves to progressively 
lower frequencies compared to twice the peak frequency 
of the tunneling data. 

One more point. For s— wave case, weak coupling cal- 
culations [|l5| have shown that for an attractive zero- 
sound coupling, Rf u u{u>) has two peaks, one at uj res and 
another near 2A. We did not find indications for a two 



peak structure in R f u u (uj) . The reason is that in our case, 
the peak in R(to) is already rather broad. Furthermore, 
we found numerically that the overall shape of Rf u u(uj) 
does not change much compared to that in R(uj), i.e., 
the dip and the linear behavior at small frequencies are 
present also in Rf u u{ui) (see Fig. ||). ^From this perspec- 
tive, the shape of the Raman intensity is mostly deter- 
mined by the self-energy corrections, while the position 
of the Raman peak is determined by the resonance in the 
Raman vertex. 

To summarize, we have demonstrated that the exper- 
imentally observed doping evolution of the B\ g Raman 
intensity in cuprates can be explained by an interaction 
with spin fluctuations. We argued that for the optimally 
doped and underdoped materials our results capture all 
salient features of the experimental data in Fig.|]: (i) the 
\fui behavior of the intensity in the normal state, (ii) a 
predominantly linear low frequency behavior of R(lu) in 
a superconductor, (iii) a reduction of the peak amplitude 
with decreasing doping and a development of a dip above 
the peak, and (iv) a progressive downturn deviation of 
the Raman peak position compared to the distance be- 
tween the peaks in the tunneling density of states. Fi- 
nally, the prediction that the Raman peak survives in the 
pseudogap regime is also consistent with the data [||. 

An issue which we didn't address in this paper is the 
experimentally observed strong discrepancy between the 
Raman data in B\ g and A\ g geometries pj| . This dis- 
crepancy (e.g., different locations of the peak frequencies) 
is still unexplained [^5| and clearly calls for more theo- 
retical work on Raman scattering. 
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